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Abstract 

Various soliton-obstruction systems have been studied from analytical perspec- 
tive. We have used collective coordinate to approach the dynamics of solitons as 
they meet a potential obstruction in a form of square barriers and holes for three 
models in (1+1) dimensions, namely: \(f) A model, deformed Sine-Gordon model, and 
a model that give rise to Q-ball solution. We have shown that our approximated 
field solution is valid enough to describe the behaviour of solitons scattering off a 
potential obstruction. 

1 Introduction 

Solitons scattering from obstructions have been studied numerically {eg [U El El HI El 
E]) and have shown interesting dynamics. As we have seen in our previous study of 
soliton scattering that each model has its on dynamical features. However, all models 
have as observed numerically mutual dynamics. In all models investigated so far solitons 
have elastic scattering on barriers and inelastic scattering on holes. Furthermore, the 
scattering on barriers have revealed that the core of solitons is not important and the 
dynamics in this case resembles the scattering of a featureless point particle. Therefore, 
results obtained from the dynamics on barriers can be worked out analytically. On the 
contrary, the scattering on holes have shown some challenging dynamics. Solitons in 
all models investigated loses some their energy which makes the analytical description 
of such behaviour more complicated. Mo, solitons have shown quantum -like behaviour 
as observed in many studies and this makes the behaviour quite hard to be understood 
within a classical theory. 

In this paper, we will try to shed some light on the dynamics of various soliton- 
obstruction systems from analytical perspective and compare the analytical results with 
the ones observed in numerical simulations. To approach the dynamics of the soliton- 
obstruction systems we will use collective coordinates, ie the parameters of the field so- 
lution. We will, for simplicity, approximate the solitonic solution by only one parameter 
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which represents the position of soliton. Hence, the soliton solution can be approximated 
as <p ( x, t) « cf) ( x — X (£)), where X (t) is the position of soliton. 

Soliton field has an infinite degrees of freedom and so a more reliable approximation 
would be if we allow the kink of the soilton to vary at the obstruction, ie tp(x, t) ~ 
ip ((3 ( t) ( x — X (t))) where (3 (t) is the parameter that describes the variation of the 
kink at the obstruction. So, Far away from the obstruction at t — > ±00, (3 — > 1. But, 
the resulting dynamics is more complicated to handle. Thus, we will consider only an 
approximated field with one parameter that describes the position of the soliton. 

In following sections we will study three different soliton-obstruction systems in (1+1) 
dimensions and show how much our approximations are valid. 

2 The Central Idea of Approaching the Dynamics of 
Soliton-Obstruction Systems 

We will Consider Lagrangian densities for which a soliton solution can be obtained. A 
Lagrangian density that describes soliton-obstruction systems is given by 

C=T-V } (1) 

where T and V are the kinetic and the potential terms respectively. The obstruction is 
introduced to the Lagrangian density via a coupling in the potential term. The obstruction 
is confined in a certain region of space and it is seen by a soliton as an external potential 
hole or a barrier. 

where V = X(x) V and A = 1 + A (x). The obstruction is localised in a finite region 
of space and are either square wells or square barriers. So, in what follows, we will write 
A (x) = Xq. Thus we can describe them using Heaviside functions which are defined by: 

Therefore the potential can be written as: 

X(x) = 1 + X [H {x+ x ) - H(x- x )], (3) 

where xq is the position of the obstruction and Ao is the parameter that describes the 
magnitude (ie height or depth) of the obstruction. A study that has been conducted in 
this regard for Sine-Gordon model in [7j, has considered the obstruction as a one point 
perturbation. In our study we will over look the dynamics that span the space for which 
the obstruction is localised. The sign of Ao determines whether the potential obstruction 
is a hole or a barrier. When Ao > 0, the obstruction is a barrier and a hole when Ao < 0. 
In the following sections we will analyse various soliton-obstruction systems that would 
shed a light on the dynamics of such systems. So, we need to calculate the Lagrangian 
which is given by 

/oo 
dxC 

/oo 
dx (T — (1 + A [H (x + x )-H(x- x )}) V (x)) 
-00 

/oo rx 
dx (T-V)-Xq dxV (x) . (4) 
-OO J — XO 
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In the following sections we will build an approximate field solutions for various soliton- 
obstruction systems and give analytical description for the dynamics of these models and 
compare our analytical description with the numerical simulation work that have been 
explained in the previous papers [H EJ [6] . 



3 Model 

The Lagrangian density for the A0 4 model is given by 

£=i^-V-A(^-l) 2 



We use the ansatz 



(p(x;X)= tank k/2 ( x - X ( £)) 



where X ( t) is the position of soliton as a function of time. 

Now, the Lagrangian density, after substituting the above results, becomes 

C = ( X 2 - 2] sech 



The Lagrangian: 



V2(x- X) 

-A [ H ( x + x ) — H ( x — x )\ sech 4 ' V% ( x — X) 



C dx 



3^2 
4 

3^2 



-An 



sech 



y/2(x- X) 



= (X 2 ~2 
= (x 2 -2 

-^-A tank (V2 ( X + x )) [ sech 2 ( X + x )) + 2 
/2 

+^-A tan/i (v 7 ^ ( X - x )) [ sech 2 (V2 ( X - x )) + 2 



The potential is given by 



(5) 



(6) 



(7) 



V{X) = ^\, 
o 



tanh 



V2(X+ xo)] ( sec/i 2 [V2 ( X + x ) 
'y/2(X- x )] ( sech 2 [V2 ( X - x 



+ 2 
+ 2 



(9) 



Fig. 1 and fig. 2 show the potential as a function of the position of soliton for An = ±1 
respectively when the obstruction is located at | xq\ < 5. 

The potential exponentially suppressed outside the obstruction's region. 
Using the Lagrange- Euler equation, 44 — -4 (§^) = 0, the e.o.m is 



3V% 



X + An [, sech 4 (V2(X+ x )) - sech 4 ( X 



0. 



(10) 



The force of the obstruction on soliton is F — MX (M is the rest mass of soliton) 
and is given by 
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Figure 1: The potential raised by a barrier, A =l 




Figure 2: The potential raised by a hole, Aq=1 
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Figure 3: The force on the soliton by a barrier, A = 1 




Figure 4: The force on the soliton by a hole, A = — 1 



F = -Ao [sech A (y/2(X + x )) - sec/i 4 (V2 ( X - x ))] , (11) 
where the rest mass M rest = 

In case of a barrier, Ao > 0, then the force is repulsive as F < 0. The case of a hole 
where Ao < the force is attractive since F > 0. Fig 3 and fig. 4 show the force exerted 
by the obstruction on a soliton for Ao = ±1 when the obstruction is inserted in the space 
region, | x \ < 5. This is in agreement with the observed behaviour. In case of a barrier 
the force for the first half of the barrier is repulsive but it becomes attractive for the 
second half. In the hole it is the other way around. 

The force as can be seen in fig. 3 and fig. 4 decays out exponentially as we move 
away from the obstruction and increases, in magnitude, exponentially at the edges of 
the obstruction. Then it dies out exponentially when the soliton is inside the hole or at 
the top of the barrier. Thus, when the soliton is far away from obstruction, X — > ±oo, 
A = 0, the force is zero and so 

X = 0. (12) 

The solution is 



X = X + u t, 



(13) 
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Figure 5: The trajectories of soliton-barrier system 
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Figure 6: The trajectories of soliton-hole system 



where Xq is the initial position of soliton. 

We have solved the equation of motion ffTOl) with initial conditions that specify the 
values of the position X (0) of the soliton and its speed using fourth order Runge Kutta 
method. Fig. 5 and fig. 6 show the trajectories of the soliton starting from an initial 
position, X (0) = —15. These trajectories describe a soliton moving with a speed of 0.5 
and interacting with a barrier of height Ao=0.25 and a hole of depth -0.25. The figures 
demonstrate the validity of our approximation. We have found such agreement between 
our analytical approach and the numerical simulations in many occasions. 

We can calculate the total energy of soliton. 



E 



dL 
dX 



X - L 



V2 
6 



Ar 



tank U/2(X + x ) 
'V2(X- x 



sech 2 ( V / 2(X+ x Q )) +2 
sech 2 (y/2(X- x )) +2 



— tank 

The total energy, when the soliton is away from the potential, A = 0, is 
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(14.) 



E 



X z + 



3V2' 



(15) 



A calculated M B observed M B 



0.25 2.12131 2.12126 
0.5 2.35701 2.35696 

0.75 2.59272 2.59262 



Table 1: The calculated and observed masses of a soliton at the top of barriers, M B 



\q calculated M B observed M B 

-0.25 1.6499 1.6499 

-0.5 1.4142 1.4142 

-0.75 1.1785 1.1785 



Table 2: The calculated and observed masses of a soliton inside holes, M H 



And so when the soliton is at rest,ie X — 0, the total energy is the mass of soliton: 

Mrest = ^=. (16) 

This is, in fact, the minimum energy of the soliton i.e E m i n = M rest . 
Therefore, 

£a i^r (17) 

In order to find mass of soliton at the top of a barrier or inside a hole with X = 
we need to position the soliton where the force of the barrier or a hole on the soliton is 
minimum. The best position is when X = since the force is zero. However, one needs to 
keep in mind that a soliton is not a point particle but an extended structure where there 
will be always an interaction between the tail of the soliton and the tail of the force on 
both ends. At X = 0, the total energy which is the mass of a soliton at the obstruction 
is given by 



M = M rest + —X tanh (v^rro) (sec/i 2 (v^xo) + 2) 

= M rest + 0.943A . (18) 

Table 1 and table 2 show the numerically and calculated masses of soliton at various 
A when x = 5 for various barrier heights and hole depths. 

Thus, the tablesl and 2 confirm the excellent agreement between the analytical and 
numerical values. 

The total energy must be conserved. So, the energy away from the obstruction must 
equal to the energy during which the soliton interacting with the obstruction. Thus, if 
a soliton is moving with a velocity u and experiences an obstruction, ie a barrier, then 
energy conservation implies 
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n ^ M rest = M, (19) 



where w& is the velocity of the soliton at the obstruction and M is the mass of the 
soliton at the top of a barrier or inside a hole. Hence, we can with a good approximation 
calculate the velocity of soliton when it is crossing the region of the obstruction 



M * 2 



L rest 



(1 - u 2 ) (20) 



To calculate the critical velocity, u c , in case of a barrier we set u = and the equation 
reduces to 



U r . 



M B - 2 



(21) 



\ \M rest , 

When the soliton moves with the critical velocity, the kinetic energy approximately 
equal to the rest mass energy at the top of the barrier. 

The equation (121]) agrees with the observed values of the critical velocity of a soliton 
moving over a barrier with an error of~ 3 percent. 

Alternatively, we found that the energy is scaled up by a factor of \[\ at the top of 
the barrier, 



M B = M rest \l X = — (22) 



3^2' 

This agrees almost with the observed values. Now, the energy conservation in the 
non-relativistic limit is 



where u c is the critical velocity. Hence, the critical velocity in the non-relativistic 
limit is 




«c*\2 VA-1 (24) 



Now, we will calculate the critical velocities for different barrier heights, i.eXo > 0, using 
the above non-relativistic limit of the critical velocity, and compare them the observed 
ones. Table 3 summarizes our calculations and observations of the critical velocities for 
different barrier heights. 

From table 3, one can see that for low barrier height and obviously low critical veloc- 
ities there is an agreement between the calculated critical velocities and the numerically 
observed ones. However, for higher barrier heights there is a disagreement and the differ- 
ences grow wider as we increase the height of the barrier. We can have a full agreement 
when we use the relativistic correction to the critical velocity. 

-^sst, = M rest \fl 

1 - V 2 = 
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Ao 


u c (non — relativistic) 


observed u c 


0.125 


0.34 


~ 0.34 


0.25 


0.49 


~ 0.45 


0.5 


0.67 


~ 0.59 


0.75 


0.80 


~ 0.65 



Table 3: The critical velocities in the non-relativistic limit vs the observed ones 



Ao 


u c (relativistic) 


observed u c 


0.125 


0.33 


~ 0.34 


0.25 


0.45 


~ 0.45 


0.5 


0.58 


~ 0.59 


0.75 


0.65 


~ 0.65 



Table 4: The critical velocities in the relativistic limit vs the observed ones 
Therefore, in the relativistic limit, the critical velocity reads 



"WttV (25) 

If we recalculate the critical velocities for the same barrier heights as in Tabel 3 we 
find an excellent agreement with numerically observed values, see table 4. 

We were, also, able to produce the critical velocity curve using our approximation for 
the model by solving (flUj) using the forth order Runge-kutta method. Figure 7 shows 
the analytical critical velocity curve for a barrier of 0.25 height. The velocity of soliton 
which produces this curve is u cr = 0.421025 which is marginally less than the numerical 
value, ie u cr ~ 0.45. The system has an infinite degree of freedom and at the critical 
velocity the time of the interaction is large therefore this would contribute to the small 
difference between the two critical velocities. However, this is a further demonstration of 
the validity of our approximation. 

4 Deformed Sine-Gordon Model 

We will analyse the scattering properties of topological solitons in a class of model which is 
the generalisation of Sine-Gordon model and which has been recently proposed by Bazeia 
etl [8]. The model depends on a positive real non-zero parameter n but we will consider 
the model for its integer values as when n = 2 the model reduces to Sine-Gordon one. 
The model is constructed by a class of potentials which is given by 

2A 2 

Vn (if) = —7 tan 2 (ip) (1 - sin n (if)) , (26) 
n z 

where the model will be approximated by the following soliton solutions 
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Figure 7: The trajectory of soliton-barrier system showing the critical velocity in case of 
a barrier of 0.25 height 



(p = sin [w] n 

_ exp\2(x-X(t))\ 
w l+exp[2(x-X(t))]- 

The parameter A ,as before, is expressed in terms of the Haveiside function 
The Lagrangian density is 

£ n = i^_ U 2 _ V n (<p). (27) 
Substituting the values of <fi, (ft' and V, the Lagrangian density becomes 

2- /-, \2 



r = 

n 



2 w n (1 - w ) / ■ 2 \ 2A 2 Wn(l- w ) 

it r ~ i J~^7 - ^ X\~" ^ 

^ ( 1 — W n\ n ( 1 — W n j 



The Lagrangian is given by 



oo 
-oo 



n Z V 7 J-oo ( 1 — u; n 



(A + 2) / ^dx. (29) 



2 V/MJ I f-, 2. 

/i J - X 1 - ti)« 



It is convenient to change variables from x to r = u; » . By applying this change of 
variables the first integral in the Lagrangian [it (1291) ) becomes 



2 9 

00 w~ (1 — w) n f 1 r(l — r n ) 



-oo (i _ 2 7o (1 - r 2 ) 

The Lagrangian can be recasted as 



dr. (30) 



L n — 

n 



-(x 2 -2) S n -\ (\ + 2)JC n , (31) 
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where 



s„= / l r(1 -: w) dr, 



/C 



/o 1 — r 2 
1 r^o r (1 — r' 

n = ~ ; 

n J-xo I — 



dr. 



-xo J- — T 

After manipulating the integrals, We find that 

\\-ln2 n=l 
Sn ~ { \ n=2. 

and all other satisfy the recurrence relation 



Sn+2 — S n + 



1 



n + 2 



/C n will be expressed in terms of w± = w± (X). 

exp (—2 ( X — xo)) 



W- = 



w+ = 



1+ exp(-2(X- x ))' 

exp(-2(X + x )) 
1 + exp (-2 (X + x ))' 



(32) 
(33) 

(34) 

(35) 

(36) 
(37) 



n — 1, .., 6. are 



' (w-- In (1 + wJ) - ( w + - In (1 + 
l 



^ + - In ( 1 + wi ) - 



™+ . -3 _ln(l + ..^ 



i 

2 4 



2 4 



^ + + wi - In ( 1 + wi ) - | ^ 



+ w+ — In ( 1 + w 



— H + — 

2 ~ 4 ~ 6 



The potential of the system is 

K(X) = A (A + 2)/C r 
The force on soliton by an obstruction is 



f)JC 

F n (X) = -\ (\o + 2)j£ 



where 



n=l 
n=2 

n=3 
n=4 

n=5 

n=6 

(38) 

(39) 
(40) 
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Figure 8: The potentials raised by a barrier of height A = 1 



dlCn 

dX 



- 2 (f^)-KGrS)) 

\ ( w- (1 - «;_) - u>+ (1 - u>+)) 



(1 



w. 



1 + wi 



1+ wi 



V l+tof/ 



«;-) i + 



™ + (i - w+ ) i + 



n=l 
n=2 



n=3 



n=4 



^wl (1 - w. 



1 + w i + — 



l+w 



1 w+ (1 - u>+) 1 + tu+ + w 



18 



iw-6 (1 - w-) 1 + + wl ) - I wl (1 - iu+) 1 + «4 + iu 



n=5 



n=6 



(41) 

The Figures 10, 11 and the figures 12, 13 show the potentials and forces that Solitons in 
this model would experience at a barrier A = ±1 respectively. The potentials and forces 
are decreasing as n increasing. The potentials are getting more deformed, ie asymmetrical, 
as n increases with the exception for the case n = 2,(ie the Sine-Gordon solution). 

Using the Euler-Lagrange equation, we can determine the equation of motion for each 
case. However, we will only consider the case where n = 2 which is the Sine-Gordon 
solution. 

The e.o.m for n = 2 is 



X-A (A + 2) («;_(! 



u>_ — ul 



(1 -w + )) = 0, 



(42) 



and the rest mass for the Sine-Gordon soliton is M Test ( 2 ) 
The e.o.m. for n = 3 is 



| (J - In 2) X - h (A + 2) wl (1 -«;_)( 1 + - 



l+w: 
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Figure 9: The forces on solitons at a barrier of height A = 1 




Figure 10: The potentials raised by a hole, Xq = — 1 



13 




Figure 11: The forces on solitons at a hole, A = — 1 
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A (A + 2)\wl(l-w+)\l+w 



1 + w%, 



0. 



and the rest mass for this soliton M r , 



(43) 



W(3) = § (I - In 2). 

When the soliton is far away from the obstruction, the equation of motion reduces to 



X = 0, X = X + ut, 
where X l is the initial position of the soliton. The total energy 



(44) 



d L n ■ 
E n = — — X — L r 

dx 



= -(x 2 + 2) S , n + Ao(A + 2)/C n 
For n — 2, the total energy is 

E 2 = \ ( X 2 + 2) + ^ (A + 2) ( w_ - w + ) 
Far away from the obstruction, the energy is simply the 

4 2 



(45) 



(46) 



(47) 



If the soliton is not moving, ie X = then the total energy corresponds to the rest 
mass energy of the soliton. Therefore, 



E 2 > 



(48) 



The rest mass energy of the soliton at the top of the barrier (ie Ao > 0) or inside a 
hole(«eAo < 0) can be calculated by fixing the soliton position over the barrier or inside 
the hole. The best choice would be when the soliton is at X = because at this position 
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A = 0.4 



n calculated M B observed M B 



1 0.908 0.908 

2 0.74 0.74 

3 0.632 0.631 

4 0.555 0.554 

5 0.497 0.496 

6 0.453 0.449 



Table 5: Solinton masses at the top of a barrier of 0.4 height 



n 


calculated M H 


observed Mh 


1 


0.417 


0.417 


2 


0.34 


0.3400 


3 


0.2904 


0.2904 


4 


0.255 


0.255 


5 


0.2016 


0.2297 


6 


0.1836 


0.2097 



Table 6: Soliton masses inside a hole of -0.4 deep 



the force on the soliton by the obstruction as can be seen from the figures is the minimum. 
In this case the functions w± = w± (0) are 



w. 



(0) 



w + (0) = 



exp(2(x )) 
1 + exp(2 (x )) ' 

exp(-2 (x )) 



1 + exp(-2 (x ))' 

Since xq = 5. this would give with a good approximation W- ~ 1 and w + 
Setting X = at the obstruction the solitons masses are given by 

M n = M rest + ^y^A (Ac + 2) 
Hence, the masses of solitons at the obstruction for n = 1, .., 6. are 



(49) 



(50) 



0. 



(51) 



(52) 



'0.614 + 0.307A (A + 2) n=l 

0.500 + 0.25A (A + 2) n=2 

0.427 + 0.2135A (A + 2) n=3 

0.375 + 0.1875A (A + 2) n=4 

0.336 + 0.168A (A + 2) n=5 

[ 0.306 + 0.153A (A + 2) n=6 

Table 5 and table 6 compare between the numerical observed masses and the calculated 
ones for a barrier of 0.4 height and a hole of -0.4 deep. 
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Figure 12: The trajectories for n = 2 soliton solution over a barrier of 0.1 height, u = 0.45 

Thus, the agreement between the numerical and calculated values for the masses at the 
top of a barrier or inside a hole is perfect. This shows the validity of our approximation. 

The critical velocities for solitons-barrier system can be calculated easily as before 
using the following equation 



For a barrier of height 0.4, the critical velocity obtained numerically for ,n = 1, ..6 
solutions is ~ 0.7 [5J. Using the above equation, the critical velocity is 0.737. The 
~ 5.3 percent difference is because the dynamics of this system is very sensitive to the 
perturbation raised by the obstruction. Unlike the A<^ 4 model, in this class of models 
solitons are excited as they meet the obstruction in their way and this would contribute 
to the masses of solitons as we have discussed that in [5]. If we subtract the excitation 
energy from the masses observed numerically as we did in [5] then there will be an excellent 
agreement between them. 

We will solve the equation of motion for n = 2 soliton solution (ie fT4*2l using forth 
order Runge-Kutta method with initial conditions that specify the speed of the soliton 
and its position. 

Fig. 14 shows a good agreement between the analytical and numerical simulation 
curves for n = 2 soliton moving with a speed of 0.45 and encountering a hole of -0.1 
depth. However, fig. 15 shows that the there is a noticeable difference between the 
analytical and the numerical curves for n = 2 soliton moving with a speed of 0.5 and 
meeting a barrier of height 0.1. We found out that as we are getting closer to the critical 
velocity of the system the difference between the analytical and the numerical simulation 
curves grow wider. In the case of a barrier of 0.1 the critical velocity is 0.425 and if we 
keep increasing the velocity the difference is diminishing. The curves in the case of the 
hole, fig. 14, the speed of soliton is much larger than the critical velocity and so the curves 
are in a good agreement. 

For a barrier of 0.4 height, the critical velocity for n = 2 solution as calculated from 
(153"j) is ~ 0.7 which agrees with numerical value. By using the fourth order Lunge-kutta 





(53) 
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Figure 13: The trajectories for n = 2 soliton solution over a barrier of 0.1 height, u = 0.45, 
a = 1.40 



25 




50 100 150 200 250 300 



Figure 14: The critical trajectory for n = 2 soliton in case of a barrier of 0.4 height, 
u c = 0.692375 



method in solving the equation of motion for n = 2 we found out that the critical curve 
can be obtained with a critical velocity, u c = 0.692375 which is in perfect agreement with 
the numerical value. Fig. 16 shows the critical trajectory in case of a barrier of 0.4 height. 

5 Q-ball-Obstruction System 

Consider the field configuration 

$ = n(t) e ie{t) f{x- X{t)), (54) 

where rj and 6 are time- dependent real moduli and f (x — X (t)) is the profile func- 
tion, which does not have any explicit dependence on t, of a Q-ball with frequency u. We 
will study the dynamics of the system without any obstruction, ie A = 1. The Lagrangian 
density in (1 + 1) dimensions of our model is 

C = I|$| 2 - I|$ ;c | 2 -2|<1)| 2 + 2|<1)| 4 - |$| 6 , (55) 
Each term in this Lagrangian density is given in the following equations 
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$|4 


= ^/ 4 


$ 6 





$ = r l e ie f+iee ie f + r 1 e ie f x X, 

^m 2 = (v 2 + v 2 e 2 ) f 2 + v 2 x 2 f x , 



®x = ve fx, 

^\$ x \ 2 = v 2 fl 



The Lagrangian is 



L — J dx C 

The Lagrangian after inserting each term of the Lagrangian density is 
L = \ (V 2 + V 2 2 ) h + \v 2 X 2 f\ - \r? I x - 2r? I 2 + 2r? J 4 - 



where 



/oo 
f n dx, n = 2,4,6 
-oo 

fx d x i 

-oo 



-oo 
roo 



/oo 
fxdx. 
-oo 

Applying the Euler-Lagrange equations 

(Pi T,\ 

0, 
0, 
= 0, 

yields the following equations of motion 

(fj - V 9 2 + 477) I 2 - V X 2 I x + V I X - 8 V 3 I A + QV 5 h = 

+ 2^0 = 
18 



dL 


d 


( dL 


drj 


~di 1 


{Or], 


dL 


d 


(dL" 


~m ~ 


It ' 


{09, 


dL 


±\ 


(dL 


dX 


dt 1 


[dX 



k+2 [l) * = °- ( 71 ) 

We can have a conserved quantity using 

d (dL\ dL 

— — r- =>• — r- = COnSt. (72) 

dt\de ) do 

from which we obtain 

?? 2 # = const. (73) 

Similarly, 
from which we obtain 



const. (74) 



?7 2 X = Const. (75) 
If f] — 0, that is when the magnitude of the Q-ball does not change with time we have 

X = X - X = ut, (76) 
where u is the velocity of the Q-ball. Also, 



The energy density 



The total energy 



Q = 6- 6 = ut. (77) 



H = + - C (78) 



E= / Hdx 

J — oo 

= 1 -(v 2 + V 2 d 2 )l2 + lv 2 X 2 I x + 1 -n*I x 

+2rj 2 I 2 - 2n 2 I 2 - 2n 4 J 4 + rf I 6 (79) 

If we assume that rj = 1 and there is no change of this amplitude as the time evolves, 
rj = then the total energy becomes 

E=^9 2 I 2 + ^X 2 I X + ^I X + 2I 2 -2I 2 -2I 4 + I 6 (80) 

If the soliton is not moving ( ieX — 0) then the total energy reduces to the rest energy 
of the soliton 

E rest = l -9 2 I 2 + l -I x + 2I 2 -2I 2 -2h+ I & (81) 

Now, let us introduce an obstruction to the system (ie Ao 7^ 0) and look over the 
dynamics of the system. Let us assume that the scattering with an obstruction that 
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leaves the Q-ball stable will not affect the magnitude of the Q-ball, ie rj (t) ~ 1. Thus, 
we can have the field configuration approximated by the phase and position parameters. 
Thus , we have 



${x- X(t)) = e ie(t) f(x- X(t)). (82) 
The Lagrangian density with an obstruction introduced to the system is 



C = ^|$| 2 - ^\® x \ 2 - A (2|$| 2 - 2|$| 4 + |$| 6 ) . 
Now, each term in this Lagrangian density is given by 

\<5>\ 2 = 6 2 f 2 + X 2 f 2 x , 



33) 



34) 



1$ I 2 = f 2 



15) 



= /» n = 2,4,6. 



The Lagrangian density becomes 

1 • 



C = §* 2 / 2 + 1 -X 2 f 2 x - 1 -f 2 x -\(2f 2 -2f+f 



The Lagrangian is 



L 



dxC 



le 2 I 2 + ^X 2 I X -^I X -2I 2 + 2I 4 - h-X (2I' 2 ~2I' 4 + I' 6 ), 



where 



r n = H dxf n , n = 2,4,6. 

J—xo 



The total energy is 



(86) 



37) 



(89) 



dX 86 
= ±X 2 I x + U 2 I 2 + ^I x + 2I 2 -2h + h 
+X (2I' 2 -2I' A + I' 6 ). 



(90) 



In case the Q-ball is far away from obstruction, the total energy reduces to AST]) . 
The field equations after using Euler-Lagrange equations for the parameters X and 9 



arc 



XIx + Xo 



ax ax ax 

= =>- = constant. 



(91) 
(92) 
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9 corresponds to the angular frequency of the Q-ball, 9 = co. 

There is an exact solution to such system in (1+1) dimension [6]. Let us approach the 
dynamics of the system by using an approximate solution given by 



f(x-X(t)) 



uj 



2 + \J2uj 2 - 4 cosh 2V4-W 2 ( x - X (t)) 



(93) 



where uj = 9. The potential of our model put restrictions on the values of cu for which 
a Q-ball solution exists. Thus , the Q-ball field exist for uj in this range, it \/2 < uj < 2. 
We will select a particular phase that simplify our model and then look over the dynamics 
of the soliton-obstruction system for this phase. We will choose, for analytical simplicity, 
uj = a/3 and the approximate solution reduces to 



f(x-X(t)) 



1 



\| 2 + V2cosh(2(x- X(t)))' 
We will evaluate the integrals,, Ix, Ix, and I n using this solution. Thus, we have 



(94) 



h = 


/ f 2 dx = 0.623, 

J — oo 


(95) 


h = 


roo 

/ fdx = 0.123, 

J — oo 


(96) 


h = 


roo 

/ fdx = 0.029, 

J — oo 


(97) 


h = 


roo 

/ fldx = 0.1885, 

' — oo 


(98) 




roo 

/ f 2 x dx = 0.1885, 

'— oo 


(99) 



Calculating the energy density away from the obstruction for a static Q-ball, (1911) . 
with u = a/3 using the above values of the integrals gives E = 2.058 which is the rest 
mass of the Q-ball and this is in agreement with the numerical value. 

Now, the Lagrangian becomes 



L = ~ (0.1885) ( X 2 - 2) - A (2 V 2 - 2 1' 4 + I' 6 ) 



(100) 



We will evaluate the integrals I' n , n = 2, 4, 6. 
For n = 2 



I' 2 = —j= tank' 1 (0.4142 tank [( X + x )}) 
V2 

— ^ tcrn/i" 1 (0.4142 tan/i [( X - x )]) 
v2 



;ioi) 



For n 
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Figure 15: The potentials raised by a barrier, A = 1 



—= tank- 1 [0.4142 tank ( X + x )} 
\J2 



-= tank' 1 [0.4142 tank ( X - x )} + 

V2 



sinh(2(X + x )) 
4V2 + 4cosh(2(X + x )) 
sink (2 ( X — xq)) 



4V2 + 4 cosh (2 ( X - x )) 



(102) 



For n = 6 



4 



htanh- 1 (0AU2tanh {X + x )) 


8cos/i (2 (X + a;o)) + Vlcosh (4 (X + x )) + 5\/2 


8 


2 + ^005/1 (4 (X + z )) 


2 



7^2smh (2 (X + xq)) - 3smft (4 (X + x )) 



2 + V2cosh(2 (X + x )) 



htanh- 1 (0AU2tanh (X - x )) 


8cosh (2 (X - x )) + v^cos/i (4 (X - x )) + 5\/2 


8 


2 + V2cosh(4: (X-xo)) 


2 



l\f2sinh (2 (X - x )) + 3sm/i (4 (X - x )) 



8 [2 + v^cos/i (2 (X-xq))] 
The potential is 

\/(X) = A (2I' 2 -2I' 4 + Q. 

Fig. 15 and fig. 16 show the potential raised by a barrier of height, A 
hole of A = — 1 depth located between | x \ < 10. 

The force between the Q-ball and the obstruction is 



(103) 



(104) 
1 and by a 
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where 



Figure 16: The potential raised by a hole, A = — 1 



F(X) 



= -A 



dV 
' ~dX 
dl' 



dI' A dl' % 



dX dX dX 



(105) 



dX 



2 + V2cosh(2(X + x )) 2 + V2cosh(2(X- x ))' 



(106) 



dX 



(2 + V2cosh(2(X+ x ))) (2 + V2cosh(2 (X - x ))) 



2 > 



(107) 



9X 



(2 + ^005/1(2 (X + X ))) (2 + ^005/1(2 (X- x )))' 



(108) 



Fig. 17 and fig. 18 show the force on the Q-ball due to barrier and a hole. 

As we have discussed in the previous chapter, the stability of the Q-ball is affected 
by the a deep potential hole while it is not affected for the case of high barriers. So, we 
considered only shallow holes. We will calculate the rest energy for a static Q-ball with 
co = \/3, ie X = 0, at the top of a barrier and inside a hole using (84) and compare 
it to the numerical values. To perform our calculation we fixed the position of soliton 
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0.5 




-0.5 1 ' ' ' ' ' ' ' 1 
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X 



Figure 17: The force on the Q-ball by a barrier, A = 1 




Figure 18: The force on the Q-ball by a hole, Aq = — 1 
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A 


calculated Eb 


observed Eb 


0.1 


2.1609 


2.1609 


0.075 


2.1352 


2.1352 


0.05 


2.1095 


2.1095 



Table 7: The calculated and observed rest energies of a Q-ball at the top of a barrier, E B 



A calculated E H observed E H 

-0.1 1.955 1.955 

-0.075 1.9808 1.9808 
-0.05 2.0066 2.0066 



Table 8: The calculated and observed rest energies of a Q-ball inside a hole, E H 

at the obstruction and a good choice as before is at X = where the interaction with 
obstruction is minimized. The rest energy at the obstruction is given by 

E = E rest + \ Q (2I' 2 -2I' 4 + Q (109) 
where I' n ( X; x ) = I' n (0; 10) for n = 2, 4, 6. 

I' 2 (0; 10) = 0.623, I' A (0; 10) = 0.123, Jg (0; 10) = 0.029 (110) 
Hence, inserting these values in f)109p we obtain 



E = E rest + 1.029A 

= 2.058 + 1.029A , (111) 

where 2.058 is the rest energy of the Q-ball with co = \f?>. Using this equation to 
calculate the rest energy of the Q-ball on the top of a barrier or inside a hole, we found 
a complete agreement between the calculated and the numerical values. 

Table 7 and table 8 show the calculated and numerical rest energies of the Q-ball on 
the top of a barrier and inside a hole respectively for various magnitude of the potential 
parameter, Ao- 

The critical velocity can be calculated using (114) 



E ^ 2 

J-'rest 



E 



2.058 



2.058 + 1.029A 



(112) 



For a barrier height of Ao = 0.01, the critical velocity, according to the above equation, 
for the Q-ball solution with uj = is ~ 0.1 which is the same as the numerical value [6]. 
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6 Conclusion 



we analysed the dynamics of some the soliton-obstruction systems that we have studied 
via numerical simulations. We were able to demonstrate that our approximated models 
were successfully able to explain most of the observed behaviour of solitons scattering 
off potential obstructions. We have approached the dynamics of such systems in (1 + 1) 
dimensions by using collective coordinates, ie parameters of the theory. In the case of 
topological solitons, A</> 4 model and deformed Sine-Gordon model of class I potentials, 
the position of soliton is the only parameter that has been used to approximate their 
models. More than one parameter can be used to approximate such models but this will 
result in further analytical complications. In the case of Q-ball in (1+1) dimensions two 
parameters were at least needed to explain the dynamics namely: the position of the Q- 
ball and the phase which is latter identified as the angular frequency of the Q-ball. Our 
analytical work was able to explain the observed rest masses in barriers and holes and to 
produce trajectories that matches the ones produced by simulations to a great extent and 
to find out the forces between the solitons and potential obstructions. The forces and the 
potentials exerted on the solitons due to the obstructions were plotted for each model. 
Further work can be done on building better approximated models which we will leave it 
for the future. 

Acknowledgement I would like to thank W.J. Zakrzewski for a helpful discussion. 

References 

[1] B.M.A.G. Piette, W.J. Zakrzewski and J. Brand, J. Phys. A 38, 10403-10412 (2005). 

[2] B.M.A.G. Piette and W.J. Zakrzewski, J. Phys. A 40, 329-346 (2007). 

[3] Kurosh Javidan, J. Phys. A 39 10565-10574 (2006). 

[4] J.H. Al-Alawi, W.J. Zakrzewski J. Phys. A 40, 11319-1131 (2007). 

[5] J.H. Al-Alawi, W.J. Zakrzewski, J. Phys. A 41, 315206 (2008). 

[6] J. Al-Alawi and W. Zakrzewski, J. Phys. A42, 245201 (2009). 

[7] Kurosh Javidan, Phys. Rev. E78, 046607 (2008). 

[8] D. Bazeia, L. Losano, J. M. C. Malbouisson, R. Menezes, Physica D237, 937-946 
(2008). 



26 



